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(54) Method and apparatus to compensate for gas chromatograph column permeability 

(57) In a gas chromatograph system having a col- 
umn with a stationary phase and a carrier gas moving 
through the column to contact the stationary phase, the 
system being useful for detecting analytes in a sample, 
a method and apparatus for predicting the retention 
times for the analytes under various conditions, the 
method utilizing the steps of: 



a) detecting the analytes under a number of sets of 
given conditions; 

b) calculating values for various parameters char- 
acteristic of the system based on -a mathematical 
model that includes the permeability of said column 
to said carrier gas; 

c) entering into the model the values for the charac- 
teristic parameters and at least one further set of 
conditions; and 

d) using the model to predict retention times for 
conditions other than those of step a; 

the apparatus composed of a data handling system 
including a model for predicting retention times at a vari- 
ety of conditions of operation of the system, the model 
having a first set of inputs including retention times 
determined during a number of operating condition 
sets; and a second input representative of the permea- 
bility of the column to the carrier gas, the data handling 
system possessing a prediction mechanism for putting 
into the model a set of conditions other than those of a) 
above and for pulling in the set of parameters, and for 
predicting retention times for those conditions other 
than those of a). 
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will I h Jnr!o Ve S° n re ' + ateS t0 935 chromat °9 ra P hic instruments and particularly to the prediction of how retention times 

2e SI Z SVS T r ra !| nQ C ° nditl0nS SUCh 85 C0 ' Umn ^mperature'carrier gas inlet aSSSSS 

sure, carrier gas flow velocity and column dimensions. ^ 

BACKGROUND OF THE INVENTION 

* ^ S n Chr °H mat09raphy inVO ' VeS ph ^ sica,, y separating constituents of a sample in a carrier gas which flows throuah 
iSSr^TT 8 COnstituente - A pulse ° f the sample is injected into i!Wli?SSjS 

n^int. S ', ntemCt W ' th 3 Stati0Pary Phase material in a column - ^ the end of the column the IndivS com 

ponents are more or less separated in time. Detection of the carrier with the separated constituents ?ol"des a oaSm 
o re tent-on times, which, by calibration or comparison with known samples, indicates EEKST^ES 
SSEST? ^r 4 " 31 "^- ^ main COm P° nente ° f ^ « system are the column, an injector wlh a m ixino "cham 
foTt !12 TJ 3 , 6 SamP ' e ' nt0 the C3rrier ' 3 deteCt0r af the outlet end of the column, controls and a computer 
Sent^L^A P ay,n9 T* ° f * he deteCt ° r - The S eneral| y shows th * ^9ht of each peak v^es iS 

im™ ^ h ° Ven 9enera " y iS USSd t0 el6Vate tem P e ^ure to maintain the sample in a volatile ZtelZ to 
sTtTooO d,SCnminat,0n ° f co ^tituents. A typical gas chromatographic system is disc.osed in U S J patent Na 

It is often desirable to be able to predict how retention times will change in response to chanaes in the column tpm 
: perature and inlet carrier gas pressure. Several cases of this are listed below 

' a nd nZ ^"f^r 006 chromato 9 ra P h y ^stem may want to utilize a method developed on another system 

and get retention tomes the same, or nearly the same, as those obtained on the other system for the same sample con 

SETS, f tW ° 935 chromat °9 raphs are ™y accurately calibrated, this nas generalfy noTbee 

due to d.fferences in column geometry from one column to another of the same type 

valuf fer n Shfl y ' 3 ^ Want *° SimU ' ate 6ffeCtS ° f Chan9es in 'mature and pressure in order to optimize their 
values for achieving satisfactory separation between the retention times of the various possible sample constlenSIn 
he mmimurn analys.s t.me. A means of achieving this is also described in the accompanying pSenUppliS I add 
ton to h,s invention, there are products on the market which also claim to be able to do Ms a PP |,catl0n - ln ««'- 

a n n lt T^Z ^ W3nting t0 be able t0 make such Pactions is to transfer a method from one column to 

another of a different geometry, possibly using a different carrier gas. and get the same retention timTpattern Means 
for achieving this are described, for example, in U.S. patent 5 405 432 

irM^i^o!^^ n eSirabl6 t0 h b l ab,e *° Pr6diCt 6ffeCt ° f Chan9eS preSSUre on retention «™ a 
invo ve changes in column outlet pressure, which is most commonly atmospheric pressure One of these situation.; 

limes' Z 5 T rrin ° ,IUCtUati ° nS atm ° SpheriC preSSure - Compensation JoSpr^u^ JZESE 
.mes more constant from one run to the next on a given system. U. S. patent 5 476 000 describes a means L !I1 
hi. Another situation is the use of the same method on two different Sc systems m fE5^^SS££2 

onThetwTs^ 

Another type of situation involving outlet pressure is the transfer of a method developed with one detector to svs 

2£ZS£2Tt de . teCt0r0perated at a diffe ™' P— ure. For examp.e, a method may be devetopedts ng a m ass 

t IZ m9 3 neai ' Z6r ° 0Utlet Pl ' essure ' and 1,16,1 used ™«™\y on systems with atmospheric outlet 

that wt, Jmnfn 1 T * * to * t0 aCCUrate * Ca,culate tha change in opTatog condSnt 

that will compensate for the change or difference in outlet pressure g conditions 

♦ J^^'T 935 h0 ' dUp timS iS thS time * takes a smaN se 9 ment °* carrier gas to flow from one end of the column 
o the other. At constant temperature, the constituent retention times are proportional to STSartto aas holio ST 
Wrth : temperature prog^m ming. there is a similar bu, more comp,ex dependence on Ve °hS£Z Became of Z' 

times accompanying changes in column temperature and inlet pressure. 

SUMMARY OF THF IMVPMTinM 

According to the teachings of this invention there is disclosed herein a method for more accurately nredictinn r*t*n 
tion times for various conditions in a gas chromatographic system utilizing the steps T 9 

a) detecting the analytes under a number of sets of given conditions; 
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b) calculating values for various parameters characteristic of the system based on a mathematical model that 
includes the permeability of said column to said carrier gas; 

c) entering into the model the values for the characteristic parameters and at least one further set of conditions; and 

d) using the model to predict retention times for conditions other than those of step a. 

There is further disclosed an apparatus composed of a data handling system including a model for predicting reten- 
tion times at a variety of conditions of operation of the system, the model having a first set of inputs including retention 
times determined during a number of operating condition sets; and a second input representative of the permeability of 
the column to the carrier gas, the data handling system possessing a prediction mechanism for putting into the model 
a set of conditions other than those of a) above and for putting in the set of parameters, and for predicting retention 
times for those conditions other than those of a). 

The accepted theory for making predictions of the retention time is based on the Poiseuille theory of flow of carrier 
gas. Through a combination of experimental measurements of carrier gas holdup times and computer simulations, we 
have deduced that the accepted theory is incomplete, especially in the very common case where the carrier gas is 
helium and the column is fabricated from fused silica. Correcting the theory is desirable for more accurate application 
in the examples and situations outlined above. Furthermore, it is important to be able to calculate the errors accompa- 
nying the use of the conventional theory. 

What the accepted theory does not take into account is the finite permeability of fused silica, or other column mate- 
rials, to helium or other carrier gasses, and the variation of this permeability with temperature. By incorporating this phe- 
nomenon into the theory, we are able to make more accurate predictions of how the carrier gas holdup times, and thus 
the retention times of the simple constituents, vary with temperature and pressure. This increased accuracy improves 
the methods for reducing retention time variations and fluctuations outlined above. 

The retention times are related to the system characteristics and operating parameters by a mathematical function 
including thermodynamic constants associated with interactions of the constituents with the stationary phase of the col- 
umn. The function relating retention times to those parameters such as temperature and pressure is based on theoret- 
ical relationships of a chromatographic system. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic drawing of a chromatographic system incorporating the invention. 

FIG. 2 is a flow chart of an embodiment of the method according to the invention. 

FIG. 3 is a longitudinal section of a portion of a chromatographic column as used in the system of FIG. 1 . 
DETAILED DESCRIPTION OF THE INVENTION 

The conventional expression for the gas holdup time, the time the carrier remains in the column, is written in differ- 
ent, but equivalent, algebraic forms by different authors. For an example, see G. Guiochon and C. L. Guillemin, "Quan- 
titive Gas Chromatography" (Elsevier, New York, 1 988), Chapter 2, equations (9) and (1 7). In U. S. patent 5,405,432 the 
holdup time is not expressed directly; instead, it is related to other quantities such as the carrier gas velocity at the col- 
umn outlet and the James -Martin correction factor for gas compressibilities. A convenient expression for the gas holdup 
time t D at constant temperature is 



In this expression, r| is the carrier gas viscosity at the column temperature. L is the column length, r is the inside 
radius of the column. p { is the column inlet pressure and p 0 is the column outlet pressure. 

Equation (1) is readily derived from Poiseuille's law. It is convenient for our purposes to express this law in a differ- 
ent form given by P. C. Carmen in Flow of Gas Through Porous Media" (Butterworths, London, 1956). p. 62, equation 




(i) 



(3.1) as 
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-r 2 dpi 
8tj dx 

Where x denotes position along the column, u x is the carrier gas velocity at x averaged across the column, and p 
is the pressure at x. In the standard theory it is assumed that the carrier is an ideal gas and that the carrier mass flow 
10 is the same at every point along the column. With these assumptions, the product u x p x is a constant independent of x 
In this case, equation (2) is readily solved for p x and u x and equation (1) for the holdup time is obtained from 



L 

t 0 = Ju; 1 dx 
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Because of the importance of the gas holdup time in the examples and situations listed above, we have tested the 
above theory by measuring the helium gas holdup times over a wide range of column dimensions (radius and length), 
inlet pressure, and temperature. The holdup time was taken to be the retention time for an air sample measured with a 
thermal conductivity detector. The results were compared with those expected from equation (1) using viscosities for 
helium as a function of temperature provided by S. J. Hawkes, Chromatographia, 37, 399, 1993. 

The measured gas holdup times are expected to vary as the viscosity when the temperature is changed at constant 
inlet pressure. Deviation from this expected behavior is observed in fused silica columns with helium as the carrier gas. 
Under most practical combinations of column geometry and operating conditions the deviations are moderate: the gas 
holdup times typically increase ten percent less than expected. However, with long columns and low carrier gas veloc- 
ities the deviations can be much greater. In this case, the holdup time increase significantly less rapidly than expected 
with increasing temperature at moderate temperatures Beyond a critical temperature, the helium flow from the outlet 
end of the column actually stops. 

When gas holdup times are measured as a function of inlet pressure at constant temperature, this phenomenon 
shows up as deviation from the pressure dependence of the Poiseuille law. As pressure is decreased from a high value, 
the gas holdup time first rises less rapidly than expected. Below a certain pressure, it rises more rapidly than expected' 

We have deduced that the deviations of our experimental results from the conventional theory are due to carrier 
gas permeating out through the column walls. A study of the literature shows that fused silica, the most common mate- 
rial for GC columns, has a high permeability to helium relative to that of other glasses. This is shown in data published 
by F. J. Norton, J. Am. Ceram. Soc. 36, 90, 1953. According to Norton the rate of carrier gas permeation through a layer 
of permeable material can be expressed as 

F = KPcmHg A cin ,/t mm (4) 



F is the rate of flow through the material in cm 3 at standard temperature and pressure (STP). K is the permeation 
velocity in cm gas (STP) mm thickness /s/cm 2 area/cm Hg (gas pressure difference). P cmHg 's the carrier gas pressure 
difference across the material in cm Hg, A cm2 is the area of the exposed material in cm 2 and t„, m is the thickness in mm 
K is strongly temperature dependent. Its value as a function of temperature can be expressed as: 



log I0 K = - 6.768 - 848.1 / T - 29900 / T 2 



(5) 
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It is convenient for our purposes to use the ideal gas law to express the loss of carrier gas in terms of a change in 
Px u x. change the pressure units to Pascals, take into account the cylindrical geometry of the column, and express 
length in meters and radius in mm. Making these changes we obtain 

d(u,p x ) = -IjSJL PsP» jx 

133 T s ln(r/r 0 ) r 2 (6) 



d (u x Px) 'S the change in u x p x in units of a Pa m/s over a length of dx m. T s is the normal temperature in K, P s is 
the normal pressure in Pa, and r 0 is the outer radius of the column. 

In order to calculate u x , equations (4) and (5) must be solved simultaneously. These equations are readily inte- 
grated numerically using a modem laboratory or office PC. The u x values thus calculated are then entered into equation 
15 (3) to calculate the gas holdup time. In this case, equation (3) cannot be solved in closed form, but must be integrated 
numerically. 

For isothermal chromatography, the modified calculation of carrier gas holdup time is the only change required to 
calculate retention times. For temperature programming, the situation is not this simple. Retention time t Rl for a given 
constituent of the sample is calculated from the integral equation: 



Ju x (l + f e b/T + eT )- ! dt = L 



(7) 



where a, b and c are thermodynamic parameters for the constituent and its interaction with the stationary phase mate- 
rial. 0, the phase ratio of the column, is equal to r/2 d, where d is the thickness of the stationary phase. The.carrier gas 
velocity u x is a function of temperature, which is in turn a function of time when temperature programming.; In solving 
30 equation (7), one must keep track of the position x of the constituent as a function of time. 

In the conventional theory, the average carrier gas velocity along the column at temperature T, u avg = L/t 0 , is sub- 
stituted for u x in equation (7). The retention times are then calculated from equation (8) : 



35 



ju avg (l + f e " T " T )-'cft = L 



40 

This is easier to solve than equation (7) because u does not depend on x and it is not necessary to keep track of 
the position of the constituent. Within the framework of the conventional theory this is valid because the pattern of the 
velocities along the column does not change with temperature. When permeation is taken into account, as described 

45 above, the pattern of the velocities does change with temperature. Because of this equation (7) must be used. 

The above does not consider the effects of the stationary phase or the column's polyimide outer coating on the per- 
meation of carrier gas through the column. There are two possible approaches to these layers. One is to assume that 
their impedance to carrier permeation is negligible compared to the fused silica and thus ignore them. The other is to 
determine what they are, either through a literature search or experiment, and to include them if necessary. 

50 Another aspect of the invention involves the use of other column materials and coatings to minimize carrier gas per- 

meability. We have experimental data on a very long stainless steel column that indicates that a small, but noticeable, 
amount of helium permeation is occurring. It is contemplated that an aluminum coating, which is sometimes placed on 
fused silica columns, would significantly impede carrier gas permeation. 

Norton (reference above) shows that many different glasses exhibit significantly less helium permeation than fused 

55 silica. It is part of this invention that other less permeable glasses be used to reduce or eliminate carrier gas permea- 
tion. One possibility is to fabricate composite columns with an inside surface of fused silica surrounded by a layer of less 
permeable glass. This would preserve the desirable inside surface properties of fused silica while significantly reducing 
carrier permeation. 
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The above invention may be utilized in several ways in the examples and situations outlined in the Background of 
the invention. The most straightforward is the use the theory directly as described. This involves using equations (2) and 
(6: solve for the u x values as a function of both x and T, and to use these values to solve equation (7) for thermody- 
nan\ - properties. These thermodynamic properties are then used to solve equation (7) again to determine standard- 
izec retention times for a desired set of conditions. This requires assuming that the nominal geometric parameters of 
the column, which is all that is generally known, are close enough to the actual for good accuracy. While imperfect, this 
will produce better accuracy than the conventional theory. 

There are also any number of approximate approaches which may be used. For example, with the procedures 
described above for correctly calculating holdup times, the holdup times for a wide range of situations can be simulated. 
The results can then be used to deduce a correction to equation (1) along with the range of its validity. This can be used 
in equation (7) to estimate retention times. While this procedure is only approximate, it has the advantage of requiring 
only minor modification to existing programs. The invention is utilized in a conventional gas chromatographic system 
such as described in the above mentioned U.S. patent No. 5,476,000. 

In Figure 1, a typical GC system 10 is shown. A carrier gas from a pressure-regulated source 12 is supplied to an 
injector device 14. A typical carrier gas in a GC system might be argon, helium, hydrogen, methane or nitrogen. A por- 
tion of the carrier is passed from the injector 14 into and through a chromatographic column 16. A pressure transducer 
18 is connected to measure pressure at the inlet to the column. 

A pulse of the sample is injected into the carrier in the injector device where a mixture is formed with the carrier 
gas. The sample used for characterization may contain organic molecules in combination with chlorine, oxygen, nitro- 
gen or sulphur. The pulsed mixture passes through the column, typically taking several minutes. A stationary phase of 
a suitable substance on the inside column wall interacts with the chemical constituents of the sample. Different constit- 
uents have different affinities for the stationary phase and thereby exit the column at different characteristic retention 
times. The velocity of the carrier gas contributes to the total retention time. 

A second pressure transducer 20 is connected to measure pressure at the outlet of the column. A detector 22 at 
the column outlet measures a physical property of the carrier and mixture, the magnitude of the property changing with 
each constituent passing through. Various types of detectors are used, such as hot wire, flame ionization, electron cap- 
ture, thermionic and flame photometric. The detector effects signals on a line 24, the signals being representative of the 
retention times. The column 16 is enclosed in an oven 26. 

A computer 28 regulates the pressure of carrier gas source 12, receives pressure information from inlet transducer 
18, controls the temperature of oven 26, receives pressure information from outlet transducer 20 and also receives and 
processes the signals on line 28 into a series of peaks representative of the sample constituents, the plotted locations 
of the components representing corresponding retention times. The inlet pressure is controlled by controlling the pres- 
sure of the carrier gas source. The carrier gas source can be controlled by a pressure program which could maintain a 
constant pressure or cause the pressure to vary. The oven can be controlled by a temperature program which could 
include an isothermal method or one that varies the temperature by controlling starting and final temperatures, ramp 
rates and dwell times at each temperature. The column outlet pressure can also be controlled by a pressure program 
which could maintain a constant pressure or cause the pressure to vary. The carrier gas source, inlet pressure and out- 
let pressure can also be controlled by a carrier gas flow velocity program which can cause a constant or variable carrier 
gas flow velocity within the column. The computer displays the retention times directly or other indicators computed 
from the retention times. The components are then identified by an operator or the computer as known chemical con- 
stituents, with the peak heights providing a quantitative measure. 

Figure 2 illustrates one embodiment of the method of the invention. Primary system 10 is a typical gas chromatograph 
as shown in Figure 1 . The first step is to supply the system with a test sample 30 and enter operating parameters includ- 
ing an inlet pressure 32 and various temperature programs 34. The system then operates 36 and generates retention 
times 38 which are then identified 40 for each constituent of the test sample under each temperature program. Outlet 
pressure 42 is also read. The next step is to calculate the carrier gas velocities 44 at each of several different distances 
(x's) along the column at each of several temperatures over the anticipated range of system operating temperatures. 
We have determined experimentally that the number of positions x, may be a number greater than 200. The result is a 
two-dimensional array of values. Solving 46 for these velocities requires the column dimensions 48, the permeability 50 
of the column material to the carrier gas as a function of pressure, the inlet pressure 32, the outlet pressure 42, as well 
as equation (2) 52 and equation (6) 54. The velocities 44 along with the temperature programs 34 and the retention time 
38 for each sample constituent 40 under each temperature program are entered into equation (7) 56. Equation (7) is 
solved for each of three different sets of data. The thermodynamic parameters, a, b and c are then calculated 58 for 
each sample constituent. The thermodynamic constants a and b are related to enthalpy and entropy and are tempera- 
ture dependent. Thermodynamic constant c is introduced to compensate for this temperature dependency. All of a, b 
and c are deemed to be constant for each sample constituent but generally are different for different constituents and 
stationary phases. The thermodynamic constants, a, b and c are then used to solve equation (7) 60 again as part of a 
model to predict retention times for other operating conditions 61 . 
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Figure 3 shows the construction of a column. A tube 62 is constructed of fused silica or other suitable material. The 
appropriate stationary phase 64 is present on the tube's inside surface. A coating 66 is present on the outside surface 
of the tube to impede permeation of the carrier gas through the tube. The coating layer is only as thick as required to 
obtain the desired permeation impedance to the carrier gas. Example coating materials are aluminum, glass or stain- 
5 less steel. 

Claims 

1 . In a gas chromatograph system having a column with a stationary phase therein, and a carrier gas moving through 
io said column to contact said stationary phase said system being useful for detecting analytes in a sample injected 

into said system, a method for predicting the retention times for said analytes under various conditions, said 
method comprising the steps of: 

a) detecting said analytes under a number of sets of given conditions; 
75 b) calculating values for various parameters characteristic of the system based on a mathematical model that 

includes the permeability of said column to said carrier gas; 

c) entering into said model said values for said characteristic parameters and at least one further set of condi- 
tions; and 

d) using said model to predict retention times for conditions other than those of step a. 

20 

2. The method of claim 1 wherein said sets of given conditions include temperature programs for said column. 

3. The method of claim 2 wherein said temperature programs are isothermal. 

25 4. The method of claim 2 wherein said temperature programs comprise starting and final temperatures, ramp rates 
and dwell times at each temperature. 

5. The method of claim 1 wherein said sets of given conditions include column inlet pressure programs. 

30 6. The method of claim 5 wherein said column inlet pressure programs maintain constant column inlet pressure. 

7. The method of claim 5 wherein said column inlet pressure programs cause column inlet pressure to vary. 

8. The method of claim 1 wherein said sets of given conditions include carrier gas flow velocity programs. 

35 

9. The method of claim 8 wherein said flow velocity programs maintain constant carrier gas flow velocity. 

10. The method of claim 8 wherein said flow velocity programs cause carrier gas flow velocity to vary. 

40 11. The method of claim 1 wherein said sets of given conditions include column outlet pressure programs. 

12. The method of claim 11 wherein said column outlet pressure programs maintain constant column outlet pressure. 

13. The method of claim 1 1 wherein said column outlet pressure programs cause column outlet pressure to vary. 

45 

14. The method of claim 1 wherein said sets of given conditions include column dimensions. 

15. The method of claim 14 wherein said column dimensions include column length, column diameter and thickness of 
said column stationary phase. 

50 

16. The method of claim 1 wherein said various parameters characteristic of the system include thermodynamic char- 
acteristics. 

17. The method of claim 16 wherein said thermodynamic characteristics include factors for entropy and enthalpy of 
55 said carrier gas and a factor to account for changes in said entropy and enthalpy related to temperature. 

18. The method of claim 1 wherein said at least one further set of conditions includes temperature programs for said 
column. 
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19. The method of claim 18 wherein said temperature programs are isothermal. 

20. The method of claim 18 wherein said temperature programs comprise starting and final temperatures, ramp rates 
and dwell times at each temperature. 

21. The method of claim 1 wherein said at least one further set of conditions include column inlet pressure programs. 

22. The method of claim 21 wherein said column iniet pressure programs maintain constant column inlet pressure. 

23. The method of claim 21 wherein said column iniet pressure programs cause column inlet pressure to vary. 

24. The method of claim 1 wherein said at least one further set of conditions include carrier gas flow velocity programs. 

25. The method of claim 24 wherein said flow velocity programs maintain constant carrier gas flow velocity. 

26. The method of claim 24 wherein said flow velocity programs cause carrier gas flow velocity to vary. 

27. The method of claim 1 wherein said at least one further set of conditions include column outlet pressure programs. 

28. The method of claim 27 wherein said column outlet pressure programs maintain constant column outlet pressure. 

29. The method of claim 27 wherein said column outlet pressure programs cause column outlet pressure to vary. 

30. The method of claim 1 wherein said at least one further set of conditions include column dimensions. 

31. The method of claim 30 wherein said column dimensions include column length, column diameter and the thick- 
ness of said column stationary phase. 

32. The method of claim 1 wherein said carrier gas is a gas selected from the group consisting of argon, helium, hydro- 
gen, methane and nitrogen. 

33. In a gas chromatograph system including a temperature programmable oven, a column positioned inside said oven 
having an inlet and an outlet, means for passing carrier gas through said column, means for controlling the pres- 
sure at said inlet, means for measuring the pressure at said outlet, injection means for injecting a pulse of sample 
having constituents into said carrier gas at said inlet to effect a mixture passing through said column subject to 
retention times for said constituents, detector means receptive of said constituents for effecting signals represent- 
ative of peak heights verses retention times of said constituents and processing means receptive of said signals for 
presenting peak heights verses retention times, a method for predicting retention times under various conditions 
comprising: 

specifying a plurality of temperature programs for said oven; 
specifying an inlet pressure; 

injecting a sample having constituents during each temperature program; 
measuring retention times of said constituents for each temperature program; 
measuring an outlet pressure; 
determining dimensions of said column; 

determining column material permeability to said carrier gas as a function of pressure; 
dividing length of said column into a plurality of equal distances; 
selecting a plurality of temperatures from said temperature programs; 
determining viscosity of said carrier gas at each selected temperature 

calculating carrier gas velocity at points corresponding to each of said plurality of equal distances and at each 
of said plurality of temperatures as a function of said inlet pressure, said outlet pressure, said column dimen- 
sions and said permeability; 

calculating thermodynamic parameters of each constituent as a function of said carrier gas velocities, said 
temperature programs, said retention times, said column dimensions and said viscosities; and 
calculating retention times as a function of said thermodynamic parameters. 

34. The method of claim 33 wherein said carrier gas is a gas selected from the group consisting of argon, helium, 
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hydrogen, methane and nitrogen. 

35. The method of claim 33 wherein said plurality of equal distances is greater than 200. 

36. The method of claim 33 wherein said thermodynamic parameters are related to entropy and enthalpy of said con- 
stituents. 

37. The method of claim 33 wherein said column comprises a stationary phase having an effective phase thickness 
and said column dimensions comprise said phase thickness, column length, and column radius. 

38. In a gas chromatograph system having a column; a stationary phase within said column and a carrier gas for car- 
rying a sample containing analytes past said stationary phase, the improvement comprising: 

a data handling system including a model for predicting retention times at a variety of conditions of operation 
of the system, said model having: 

a) first inputs including retention times determined during a number of sets of conditions of operation of 
said system; and 

b) a second input representative of the permeability of the column to said carrier gas, 

said data handling system having prediction means for inputting to said model a set of conditions other than 
those of a) and for inputting said set of parameters, and for predicting retention times for said conditions other 
than those of a). 

39. The method of claim 38 wherein said numbers of sets of conditions of operation include temperature programs for 
said column. 

40. The method of claim 38 wherein said numbers of sets of conditions of operation include column inlet pressure pro- 
grams. 

41. The method of claim 38 wherein said numbers of sets of conditions of operation include carrier gas flow velocity 
programs. 

42. The method of claim 38 wherein said numbers of sets of conditions of operation include column outlet pressure pro- 
grams. 

43. The method of claim 38 wherein said numbers of sets of conditions of operation include column dimensions. 

44. The method of claim 38 wherein said conditions other than those of a) includes temperature programs for said col- 
umn. 

45. The method of claim 38 wherein said conditions other than those of a) include column inlet pressure programs. 

46. The method of claim 38 wherein said conditions other than those of a) include carrier gas flow velocity programs. 

47. The method of claim 38 wherein said conditions other than those of a) include column outlet pressure programs. 

48. The method of claim 38 wherein said conditions other than those of a) include column dimensions. 
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